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We consider the evolution of perturbed cosmological spacetime with multiple scalar fields in Einstein 
gravity. A complete set of scalar-type perturbation equations is presented in a gauge-ready form, and 
we derived the closed set of second-order differential equations in several useful forms. Conserved 
behaviors of the perturbed three-space curvature in the comoving gauge, ip v , under several conditions 
are clarified. Under the slow-roll conditions, the adiabatic and isocurvature modes decouple from 
each other, and in the large-scale limit we have (i) the adiabatic mode is generally conserved, (ii) for 
a couple of special potential the isocurvature modes decouple from each other and are described by 
conserved quantities, (iii) in the two field system, the isocurvature mode is described by a conserved 
quantity for the general potential. 
PACS numbers: 98.80.Cq, 98.80.Hw, 98.70.Vc, 04.62+v 



I. INTRODUCTION 

The cosmological perturbation theory [Q , handling the 
evolution of linearly perturbed spacetime and matter 
in otherwise spatially homogeneous and isotropic world 
model, plays central role in the modern theory of large- 
scale structure formation in the universe. The linear- 
ity assumption of the structures seems to hold in the 
early evolution stage and in the large scale evolution of 
our presently observable patch of the universe. Besides 
the quantitative constraint on the cosmological model 
implying the flat universe, the recent observations of 
the CMBR anisotropics in the small-angular scale by 
Boomerang and Maxima-1 experiments [pf, in our opin- 
ion, confirm dramatically the validity of the basic as- 
sumptions used in the cosmological perturbation theory, 
i.e., the linearity of the relevant cosmic structures in such 
scales. 

The theoretical introduction of an acceleration phase 
in the early universe || enables us to draw a coherent 
picture of the structure generation process. The ever- 
present microscopic vacuum quantum fluctuations inside 
a causal domain can rapidly expand to become macro- 
scopic during this acceleration phase, which in turn, can 
evolve into large-scale structures in the universe. The ac- 
celeration phase is often modelled by using a scalar field 
with a relevant potential term. Thus, while the back- 
ground scalar field drives the early universe accelerat- 
ing, the quantum fluctuation of the scalar field accompa- 
nied with simultaneously excited spacetime metric fluc- 
tuations provides seeds for the large-scale structure and 
the gravitational wave. Recent development shows the 
possibility that observational constraints can be trans- 
lated into the constraints on, or the reconstructions of, 
the form of field potential driving the inflation. Such a 
reconstruction program is popular in the context of an 
inflation driven by a single minimally coupled scalar field 
with the slow-roll conditions B . 



The early universe scenarios with multiple episodes of 
inflation have recently attracted much attention in the 
literature. Such scenarios apparently allow more freedom 
in designing the spectra of the large-scale cosmic struc- 
tures in their generation and evolution stages M . For the 
benefit of future studies, in this paper we would like to 
present the perturbed scalar field equations valid for an 
arbitrary number of scalar fields in the gauge-ready form, 
derive closed form equations in several useful forms, and 
derive general solutions in the large-scale limit. We will 
find that in various generic situations the evolution of 
perturbations can be characterized by certain conserved 
quantities: the adiabatic perturbation characterized by 
the perturbed curvature variable in the comoving gauge, 
ip v , is conserved in various situations, and the isocur- 
vature perturbations defined by a relative perturbation 
among perturbed fields, <5 </>/,/, are often described by 
some conserved quantities, as well. We start by intro- 
ducing our notations, general equations, and strategy in 
M We set c = 1. 



II. PERTURBED WORLD MODEL AND BASIC 
EQUATIONS 

We consider a gravity with an arbitrary number of 
scalar fields minimally coupled with the gravity. As the 
Lagrangian we consider 
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R is the scalar curvature, and a, 6, . . . are spacetime in- 
dices. (j) 1 is the 7-th component of N scalar fields with 
I, J, . . . — 1,2, ... ,N; I, J indices are raised and low- 
ered by Sjj. V((f> K ) is a general algebraic function of the 
scalar fields, i.e., V(<j) ) = V(cj) 1 , . . . , cj> ). Variations 
with respect to g a b and 4> lead to the gravitational field 
equation and the equation of motion: 
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where Vj = dV/(d(j> 1 ). 

We consider the general perturbations in the FLRW 
world model. In the spatially homogeneous and isotropic 
background we are considering, to the linear order the 
three different types of perturbation decouple from each 
other and evolve independently. The presence of multi- 
ple number of minimally coupled scalar fields does not 
affect the graviational wave of the perturbed spacetime 
and the the rotational perturbations of the additional flu- 
ids, directly. Thus, in the following we will consider the 
scalar-type perturbation only. As the metric we take 



ds 2 = -a 2 (1 + 2a) drf - a 2 /3 :a dndx a 

dx a dx fi , 



.,(3) 



9*0 (1 + M + 27,«|/9 



(4) 



where a(rj) is the cosmic scale factor, a, j3, 7 and p 
are spacetime dependent perturbed order variables. A 
vertical bar | indicates a covariant derivative based on 

(3) 

g a L Considering the FLRW background and the scalar- 
type perturbation, the energy-momentum tensor can be 
decomposed into the effective fluid quantities as 



T$ = (p + 5p)5<$ + -K%. 
To the background order we have: 
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where an overdot indicates a time derivative based on t. 
To the perturbed order we have: 
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The equations for the background are: 
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where K is the sign of the spatial curvature, and H = 
a/ a. The first two equations follow from eq. (|2|), and the 
last one follows from eq. (||). 

Now, we present a complete set of equations describing 
the scalar-type perturbation without fixing the temporal 
gauge condition, i.e., in the gauge-ready form: 
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(v a V,i-\s%£j{x + Hx-a-<p)=rf, (12) 
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where x = a(/3 + a7) and k is the comoving wavenumber; 
V Q and A are the covariant derivative and the Lapla- 
cian based on g^L Equations (P-^4[) are: the defini- 
tion of k, ADM energy constraint (G[j component of the 
field equation), momentum constraint (G° component), 
ADM propagation (G^ — ^SpG^ component), Raychaud- 
huri equation (G7, — G[J component), and the scalar fields 
equations of motion, respectively. 

In the following we briefly summarize the gauge-ready 
strategy suggested and elaborated in [||,f?J. Under the 
gauge transformation x a = x a +£ a with £ = a£° (0 = rj) 
the perturbed metric and fluid quantities change as [Q: 



a = «-£', (p = tp-Hi\ k = K+[m+^)S, t 1 
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As the temporal gauge fixing condition we can impose 
one condition in any of these temporally gauge depen- 
dent variables: a = (synchronous gauge), p> = 
(uniform-curvature gauge), n = (uniform-expansion 
gauge), x — (zero-shear gauge), v/k = (comoving 
gauge), 8(f) 1 = (uniform-c/) 7 gauge), etc. By examining 
these we notice that, except for the synchronous gauge 
condition (which fixes a = 0), each of the gauge condi- 
tions fixes the temporal gauge mode completely. Thus, a 
variable in such a gauge condition uniquely corresponds 
to a gauge-invariant combination which combines the 
variable concerned and the variable used in the gauge 
condition. One can recognize the following combinations 
are gauge-invariant. 
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etc. The gauge- invariant combination 5(f) 1 , for example, 
is equivalent to Scj} 1 in the uniform-curvature gauge which 
takes (p = as the gauge condition, etc. In this way, 
we can systematically construct various gauge-invariant 
combinations for a given variable. Since we can make 
many gauge-invariant combinations even for a given vari- 
able, this way of writing the gauge-invariant combination 
will turn out to be convenient practically. Generally, we 
do not know the suitable gauge condition a priori. The 
proposal made in [|]J7| is that we write the set of equa- 
tion without fixing the (temporal) gauge condition and 
arrange the equations so that we can implement easily 
various fundamental gauge conditions: eqs. (|[]lj) are ar- 
ranged accordingly. We termed it a gauge-ready method. 



III. MULTIPLE SCALAR FIELDS 



From eqs. 
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and from eqs. (|9|, |ll"| , [jl| ) we have tp v = (H 2 /H)a v . Thus, 
using eq. (O) we can show 
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The uniform-curvature gauge is convenient for han- 
dling 8(f) 1 . Equivalently, we can use the corresponding 
gauge-invariant combination Scj) 1 ^ introduced in eq. (16). 
Assuming a flat (K = 0) background, from eqs. (|14],|9 |il|) 
we can derive 
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This equation was first derived in [|S|. Compared with 
Sep 1 equation in the uniform-curvature gauge in eq. (EQ) 
the ones in the other gauge conditions are more compli- 
cated pi. The terms inside parenthesis of the RHS of 
this equation can be written as 
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Thus, the scalar field fluctuations are generally coupled 
to each other as long as we have V ^ 0. 

We introduce following gauge-invariant combinations 

Hu = T~J F ' (19) 

which represent the isocurvature modes. If we decompose 
the effective fluid quantities in eqs. (^0) in terms of the 
sum over individual ones as S[i = ^2jSfij and \x + p = 
Yljil^J + Pj) we can show Su = a 3 (S(f>u/a 3 )' where 
Sij = 6ni/(lJLi +pi) — S/jtj/Juj + pj) is the well known 
isocurvature perturbation JlOfl . From eq. (u/m we can 
derive 



-Hij. (20) 
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which follows from eqs. (Q-|2|) with eq. (0). The gauge- 
invariant combination <p v represents the adiabatic mode. 

Notice that the last terms in eqs. (g2|,g3|) vanish for 
the following situations. 

(i) A single component case. 

(ii) Another trivial case is the one with V — constant. 

(iii) V.i oc cj)i for all components. This applies if we 
have 3i70 J + V' 1 = which implies ifi 1 — in eq. (g); 
meanwhile the slow-roll conditions assume (p 1 <C ZHcf) 1 . 
Thus, in the large-scale limit, ignoring the Laplacian 
term, we have a general integral form solution 



<p v (k,t) =C(k)-D(k) 
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where C and D are integration constants indicating the 
coefficients of relatively growing and decaying solutions, 
respectively. 

(iv) 5<fiij = 0. This case is more restrictive compared 
with (iii); the RHS of eq. (|20| ) vanishes as well. The 
solution in eq. ( |25| ) remains valid, and additionally, from 
eqs. (£l|l9) we can show ipg^i = —HS^/ifi 1 



<Pv 



Thus, under the conditions (i)-(iv) ip v is conserved in 
the large-scale limit even in the multi-component situ- 
ation, see eq. (p5|). We note that the gauge- invariant 
combination ip v which shows fundamental importance in 
handling adiabatic perturbation was first introduced by 
Lukash in 1980 Q 

In order to make eq. (P3) complete we need equations 
for 8(f>ij. From eq. (HnfJO) we can derive 
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where we used eq. (E2n in the last step. Notice that only 
tp v contributes to the isocurvature modes. Equations 
( p3| , p6| ) provide a complete set of equations in terms of the 
adiabatic, <p v , and isocurvature, 6<f>u, perturbation vari- 
ables. For a system of N scalar fields eq. fll7|) provides a 
coupled N set of second order differential equations for 
6(f) 1 . Equations (E3J,E6J) provide an alternative expression 
in terms of tp v and 8<j>u. 

In the case of two minimally coupled scalar fields, using 
ip v and 54>i2, eqs. (p3|,p6|) lead to 
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n 4>2 

Notice that the adiabatic and isocurvature modes are 
coupled generally even in the large-scale limit. Recent 
study of the two-field system can be found in [[L2| . 

IV. SLOW-ROLL SITUATION 

We have shown in the paragraph surrounding eq. ( Ba ) 
that under the slow-roll condition, i.e., for Vj = —3H<f>j 
thus (pi = 0, the adiabatic mode decouples from the 
isocurvature modes; we emphasize again that this con- 
dition differs from the ordinary slow-roll condition Vj ~ 
—3H<j>i used in the inflation models. Under such a condi- 
tion the RHS of eq. (E6h vanishes, thus the isocurvature 



modes decouple from the adiabatic one, as well. Notice, 



however, that equations for 84> v in eq 
erally coupled. In such a case eqs. (P3| 
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In the large-scale limit, we have the solution for an adi- 
abatic mode 



H 
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which is consistent with eq. (pq) because J2j &' v-J — 
constant in the slow-roll limit. Therefore, under the slow- 
roll conditions, the adiabatic and isocurvature modes de- 
couple from each other, and ip v is temporally conserved 
in the large-scale limit for an arbitrary potential. 

We can derive the solution of eq. (|3y) for the following 
cases. 

(i) V((j)\ . . . ,(f> N ) = Y,k V k^ K ) : In such a case ' the 
term with summation in eq. (pQ) vanishes, and the equa- 
tion can be written as 
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thus in the large-scale limit we have a general solution 
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where Cu and Du are coefficients of the relatively grow- 
ing and decaying solutions. These behaviors are consis- 
tent with the results in jl|] where the nontransient so- 
lutions in the zero-shear gauge were derived based on 
some additional conditions; for translation between the 
two gauge conditions we can use 
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which follows from eq. (|21l) , and the gauge- invariance of 

(ii) V{4> 1 , . . . , 4> N ) = HkVk(4> k )- In this case, we can 
show that the third term with summation in eq. (|30| ) 
becomes 6H5(f>u, thus the equation can be written as 
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where we ignored quadratic order terms in H j H 2 . In the 
large-scale limit we have a general solution 
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Therefore, when we have potential of additive or multi- 
plicative forms, as in (i) or (ii) above, the isocurvature 
modes decouple from each other, and we have 5<pu/H 
or HS(f>ij temporally conserved in the large-scale limit, 
respectively. 

(iii) With an arbitrary potential, to the linear order 
in H/H 2 , the isocurvature modes are generally coupled 
to each other. However, since Vij ~ 0(H), if we ig- 
nore linear order terms in H/H 2 , eq. ( p0[ ) becomes 
a~ 3 (a 3 5<j>u)' + (k 2 /a 2 )S(j)ij = 0, and in the large-scale 
limit we have 84>u = Cjj — Djj L a~ 3 dt, thus d(f>ij's 
are conserved. 

(iv) In the two component system eq. (pq) reduces to 
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Considering that V u ~ O(H) we have V 12 -^- ~ 
O(H). Thus, ignoring quadratic order terms in H/H 2 , 



eq. (p7|) can be written in the following form 
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In the large-scale limit we have 
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Therefore, ignoring the transient solutions we have tp v 
and Scj)i2/Q remain constant, thus 6<j)i2 oc Q. We can 
show that (i) for V(<fri,<j>2) — V\{4>i) + ^2(^2) we have 
Q = H, and (ii) for U(0i,(/> 2 ) — Ui(</>i)U 2 (</> 2 ) we have 
Q = 1/ ' H . The corresponding solutions in the zero-shear 
gauge have been investigated in |u|. In [J13J by impos- 
ing additional conditions on perturbed order variables 
nontransient parts of the solutions for 54>L and ip x were 
derived. The translation into our gauge condition can 
be made using eq. (|3J), and results in [jl3[ are consis- 
tent with ours which are represented by the conserved 
behaviors of ip v and Scp^/Q. 

Aspects of the large scale structure formation with 
multiple episodes of inflation based on the multiple num- 
ber of minimally coupled scalar field have been studied 
in ||1,1,@. 



V. DISCUSSONS 

In this paper we have investigated aspects of scalar- 
type cosmological perturbation in the context of multiple 



number of scalar fields in Einstein gravity. New results 
found in this paper are the followings. We have presented 
a complete set of scalar-type perturbation equation in a 
gauge-ready form: these are eqs. (|9T4). The conserved 
behaviors of ip v under certain conditions are the central 
theme studied in this paper: see, in particular, the para- 
graph below eq. (|24|), and eqs. ( |25| , |3l| ). A set of closed 
form differential equations in terms of the adiabatic and 
the isocurvature perturbation variables, ip v and 6<pu, is 
derived in eqs. (|23j,|26j). The equations in two-component 
situation are presented in eqs. ( p7| , p8[ ). Under the slow- 
roll conditions the adiabatic and the isocurvature modes 
decouple from each other, see eqs. (^9|,[}0|), and various 
general solutions are presented below eq. (p0|). In the 
large-scale limit we have (i) the adiabatic mode is gener- 
ally conserved, eq. (31), (ii) for a couple of special po- 
tential the isocurvature modes decouple from each other 
and are described by conserved quantities, eq. (B3lpq), 
(iii) in the two field system, the isocurvature mode is de- 
scribed by a conserved quantity for the general potential, 
eq. ©. 

A complete set of perturbation equations in the case 
of generalized gravity with multiple minimally coupled 
scalar fields is presented in Sec. 4.2 of J7p in the gauge 
ready form. Our complete set of equations in the gauge- 
ready form and other closed form equations and the large- 
scale solutions derived in this paper will be useful in han- 
dling situations involving multiple number of scalar fields 
in diverse cosmological situations. 
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